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Abstract 

We performed NLO QCD fit of combined SLAC/BCDMS/NMC DIS data at high 
^^ • X. Model independent x-shape of high twist contribution to structure function F^ is 

extracted. Twist-4 contribution is found to be in fair agreement with the predictions of 
infrared renormalon model. Twist-6 contribution exhibits weak tendency to negative 
r^ ' values, although in the whole is compatible with zero within errors. 

D ■ 

'■^ ' 1. It is well known that on the basis of the operator product expansion the deep inelastic 

scattering (DIS) cross sections can be split into the leading twist (LT) and the higher twist 
(HT) contributions. To the moment the LT contribution is fairly well understood both 
from the theoretical and experimental points of view. The HT contribution is not so well 

T^ ! explored. The theoretical investigations of HT meet with the difficulties because in the 

^ \ region where HT are most important the perturbative QCD calculations cannot be applied. 

Q I There are only semi- qualitative phenomenological models for HT calculation in literature. 

^ ■ These models are often based on the phenomenological considerations and contain adjusted 
parameters, which are to be determined from the experimental data. Unfortunately, relevant 
experimental data, especially on the longitudinal structure function F^, are sparse, come from 
. different experiments and therefore are difficult for interpretation. There are estimations of 
(j_) twist-4 contribution to the structure function F2 [|l], 0, ^ obtained from the combined fit 

-fn to SLAC/BCDMS and SLAC/BCDMS/NMC data i, |, |. Extraction of HT contribution 






to -F3 was obtained in p, ^ from the fit to CCFR data [T^. These estimations are model 



^ . independent, i.e. do not imply any x-dependence of HT, and hence a phenomenological 

H \ formula can be easily fitted to them. As to the experimental data on HT contribution to 

■ " " ' Fj,, to the moment they are available only from the QCD motivated fits to world data on 

the structure function R = a^/aT- The first fit of this kind was presented in ||TT|. This 



fit was recently renewed in ||T3[ with the inclusion of new data from the experiments E-143 



and E-140X. The world data on R were also analyzed using a QCD based model with the 



account of HT contribution [|IJ]. As far some models predict the HT contribution to the 
structure function F^, the comparison of these models with the data requires extraction of 
the HT contribution to F^ from the data on the HT contribution to R and F2. This causes 
the problems with interpolation between data points and error propagation. In addition, the 
HT contributions to R, obtained in ah these fits are model dependent, i.e. a priori suppose 
a certain x-dependence of HT. 

2. In this paper we present the results of the analysis of DIS data aiming to obtain the 
estimation of the model independent HT contribution to F^. The work is the continuation 
of our previous analysis |^, where estimation of the HT contribution to F2 was obtained. 
Our consideration is limited by the region of a; > 0.3, where the non-singlet approximation 
is valid. A data cut x < 0.75 was made to leave the region where deuteron effects are small. 
At the beginning the ansatz used in this work is essentially the same as in [0. We fitted the 



data within NLO QCD with the inclusion of target mass correction (TMC) |T5[ and twist-4 
contribution in factorized form: 

:,{P,D),HT , ^N _ -rp{P,D),TMC f ^^ \^ , "-2 \^) 



Fr^^'-^(x,g) = Fr^^'^^"-(x,g) 



1 + 



Q' 






^3/2 ^2 Q2 ^2 y^ ^2 

where F2 ' (x, Q) are the LT terms, 

i=^—^^ r = l + 



1 + v^' Q2 ' 

M is nucleon mass, x and Q^ are regular lepton scattering variables. The LT structure 
functions of proton and neutron were parametrized at the initial value of Qg — ^GeV^ as 
follows: 



Fi{x,Qo)=A,x''-{l-xf^^ 



2 

p 



F^{x,Qo)=A,,x''-{l-xf-^. 



Here conventional normalization factors Np and A^„ are 



Nr,r, = I (ixx'''''"~Vl -X) 



10 



These distributions were evolved in NLO QCD approximation within MS factorization 
scheme. The functions /i2 ' {^) were parametrized in the model independent way: their 
values at a; = 0.3,0.4,0,5,0.6,0.7,0.8 were fitted, between these points the functions were 
linearly interpolated. 

Comparing to the work we added to the analysis the NMC data (30 points on 
proton and 30 points on deuterium targets) The number of data for each experiment and 
target are presented in Table L We accounted for point-to-point correlations of data due 
to systematic errors analogously to our previous papers [^ 0, ||. Systematic errors were 
convoluted into correlation matrix 

dj = SijCXiaj + fifjisf ■ sf ), 

where vectors sf contain systematic errors, index K runs through data subsets which are 
uncorrelated between each other, i and j run through data points of these data subsets. 
Minimized functional has the form 

K,i,j 

where Eij is the inverse of Cij. Dimension of sf differs for various data set, the concrete 
numbers for each experiment are present in Table I. The factors ^k were introduced to allow 
for renormalization of some data sets. In our analysis these factors were released for the old 
SLAC experiments in view of possible normalization errors discussed in W^\. As for E-140, 



Table 1: The number of data points (NDP) and the number of independent systematic errors 
(NSE) for the analysed data sets. 

Experiment NDP (proton) NDP (deuterium) NSE 

BCDMS 223 162 9 

E-49A 47 47 3 

E-49B 109 102 3 

E-61 6 6 3 

E-87 90 90 3 

E-89A 66 59 3 

E-89B 70 59 3 

E-139 - 16 3 

E-140 - 31 4 

NMC 30 30 13 

TOTAL 641 602 47 



BCDMS and NMC data subsets, we fixed these factors at 1 and accounted for their published 
normalization errors in the general correlation matrix. 

At the first stage of the analysis we reduced all F2 data including BCDMS and NMC to 
the common value of R [|11[|. The results of the fit within this ansatz are presented in Table 



II (column 2). For the comparison at column 1 we also presented the results of analogous fit 
from from |]^, performed without the NMC data. Addition of the NMC data increased the 
value of as{Mz)-, but within one standard deviation. In the whole the figures from column 
2 are compatible with the results of earlier analysis [Q . 

The next step of our analysis was to change the form of HT contribution from factorized 
to an additive one: 

Ff •-)'-^(x, Q) = Ff '-)'™-(x, Q) + Hr\x)'-^ 

with H2 ' (x) parametrized in the model independent form analogously to /12 ' (x). We 
preferred to switch to this form because for the factorized parametrization HT term con- 
tains latent log factor and twist-6 contribution originating from F2^{x,Q) and target mass 
corrections respectively. In addition, this form is more convenient for comparison with some 
model predictions. The results of the fit with additive HT parametrization are presented 
in column 1 of Table HI and Fig.l. For the comparison we also pictured in Fig.l the value 
of F^'^'^^iXjQ) ■ h2{x) for the fit from column 2 of Table II with factorized form of HT, 
calculated at Q^ = 2.5 GeV^. One can see that the switching of the form leads to small 
decrease of HT terms at high x. Besides, their errors became smaller, meanwhile errors of 
the other parameters increased. We connect this effect with that the additive HT form is 
not so constrained as the factorized one. The change of the value of HT contribution is 
amplified due to large correlations of the HT terms with other parameters, especially with 
as (see discussion below). 

To extract F^ we replaced data on F2 by the the data on differential cross sections and 



Table 2: The results of the fits with factorized parametrization of HT. The parameters ^ 

describe the renormalization of the old SLAC data, ^2 (3 4 5 6 7 8) ^^^ ^^^ fitted values of the 

HT contribution at a; = 0.3, 0.4, 0.5, 0.6, 0.7, 0.8. For the description of the columns see text. 

1 2 

Ap 0.516 ±0.022 0.514 ±0.021 

ap 0.765 ± 0.028 0.766 ± 0.028 

bp 3.692 ± 0.032 3.690 ± 0.032 

A„ 4.8 ±4.1 4.8 ±3.9 

an 0.118 ±0.097 0.119 ±0.095 

bn 3.51 ±0.11 3.51 ±0.11 

asiMz) 0.1180 ±0.0017 0.1187 ± 0.0016 

/i^3 -0.120 ±0.017 -0.122 ±0.017 

h(^ -0.046 ±0.025 -0.054 ± 0.025 

/i^5 0.059 ± 0.043 0.043 ± 0.042 

h^Q 0.392 ± 0.076 0.363 ± 0.074 

/i^7 0.82 ±0.13 0.77 ±0.12 

his 1.54 ±0.25 1.47 ±0.24 

h^s -0.123 ±0.018 -0.125 ±0.018 

/i^4 -0.003 ±0.026 -0.012 ±0.025 

/i^5 0.162 ±0.043 0.145 ±0.042 

h^Q 0.439 ± 0.076 0.410 ± 0.073 

/i^7 0.79 ±0.12 0.75 ±0.11 

h^s 1.87 ±0.26 1.81 ±0.25 

Cp,49A 1.016 ±0.018 1.016 ±0.016 

^DA9A 1.006 ±0.017 1.008 ±0.015 

^P,49B 1.028 ±0.018 1.028 ±0.015 

^D,49B 1.012 ±0.017 1.013 ±0.015 

Cp,6i 1.021 ±0.021 1.021 ±0.019 

^Dfii 1.004 ±0.019 1.006 ±0.018 

^p'sT 1.025 ±0.017 1.025 ±0.015 

^D,87 1.012 ±0.017 1.013 ±0.015 

^P,S9A 1.028 ±0.021 1.029 ±0.019 

^D,89A 1.004 ±0.021 1.006 ±0.019 

Cp,89i? 1.022 ±0.017 1.021 ±0.015 

Cd,89B 1.007 ±0.017 1.008 ±0.015 

^D,i39 1.009 ±0.017 1.010 ±0.015 

X^/NDP 1178.9/1183 1258.4/1243 
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Figure 1: The results of fits with different forms of HT contribution (full circles - additive, 
empty circles - factorized). For the better view the points are shifted to left /right along 
X-axis. 



fitted them using the formula 



dxdy 



A-KO^ 



M2 



Q' 



{Mxyf 



F^i 



x,Q)+ 



m. 



+ (l-2^)-(F. 



y^ 



:^HT, 



X,Q) 



:^HT, 



X,Q) 



where s is total c.m.s. energy, mi is scattered lepton mass and y - lepton scattering vari- 
able. With the test purposes we first performed the fit with Fl form motivated by -R1990 



parametrization 1 1 



■,{P,D),HT 



X,Q) 



^iP,D),HT 



X,Q) 



1- 



1 + 4MV/g2 



R{x,Q) 



hi 
21n(g/0.2) 



12Q^ 



0.1252 



0.1252 + ^2 



l + R{x,Q) 



, 1 GeV^ , 
+ h2 ^, + 63 



Q' 



1 GeV^ 
Q4 + 0.32 



g2 + i 

with fitted parameters &i,2,3 and Q measured in GeV . This parametrization is constructed 
from the term with log Q-behaviour, which mimic LT contribution and HT terms with power 
Q-behaviour. The resulting values of the parameters hi = 0.100 ± 0.042, 63 = 0.46 ± 0.11, 
63 = —0.14 ± 0.18 are in agreement with the values obtained in |ll]] (61 = 0.064, 63 = 0.57, 
63 = —0.35). In our fit twist-6 contribution to R{x,Q) is compatible with zero. We can 
note in this connection, that the correlations between the parameters is large (see Table IV) 
and as a consequence the data used in the fit have limited potential both in separation of 
twist-4/twist-6 and log/power contributions to R. 

To achieve more precise determination of HT contribution to Fl one can substitute for 



the LT contribution a perturbative QCD formula |T6[ instead of phenomenological log term. 
In the leading order on a^ and for nonsinglet case it looks like 



■,{P,D),LT 



X,Q) 



27r 3^ 



1 dz 
—R 



{P,D),LT. 



z,Q). 



Table 3: The results of the fits with additive parametrization of HT. Figures in parenthesis 

are global correlation coefficients. The parameters ^ describe the renormalization of the old 

SLAC data, -f^2 (345678) ^^^ -^l'(3 45 6 78) ^^^ ^^^ fitted values of the HT contribution at 
X = 0.3, 0.4, 0.5, 0.6, 0.7, 0.8. For the description of the columns see text. 

1 2 3 

Ap 0.478 ±0.021 0.485 ±0.023 0.486 ± 0.022 

ap 0.830 ±0.034 0.816 ± 0.037 0.817 ± 0.035 

bp 3.798 ±0.043 3.791 ± 0.047 3.792 ± 0.045 

An 4.8 ±4.5 4.9 ±4.8 4.7 ±4.7 

an 0.12 ±0.11 0.12 ±0.13 0.12 ±0.12 

bn 3.58 ±0.14 3.59 ±0.15 3.58 ±0.14 

as{Mz) 0.1190 ±0.0021 0.1170 ± 0.0021 0.1170 ± 0.0021 

H^3 -0.0374 ± 0.0067 -0.0303 ± 0.0068 -0.0273 ± 0.0067 

Hi,i -0.0159 ±0.0085 -0.0031 ± 0.0084 -0.0059 ± 0.0082 

H^^ 0.0080 ± 0.0099 0.0176 ± 0.0096 0.0189 ± 0.0095 

H^Q 0.0416 ± 0.0096 0.0497 ± 0.0095 0.0494 ± 0.0092 

H^^ 0.0404 ± 0.0073 0.0529 ± 0.0077 0.0501 ± 0.0073 

His 0.0323 ± 0.0057 0.0376 ± 0.0081 0.0400 ± 0.0070 

H^^ -0.0348 ± 0.0068 -0.0205 ± 0.0066 -0.0221 ± 0.0065 

H^^^ -0.0029 ± 0.0071 0.0057 ± 0.0069 0.0067 ± 0.0068 

H^^ 0.0195 ± 0.0074 0.0274 ± 0.0071 0.0269 ± 0.0070 

H^Q 0.0330 ± 0.0069 0.0381 ± 0.0068 0.0380 ± 0.0067 

H^y 0.0301 ± 0.0050 0.0363 ± 0.0051 0.0372 ± 0.0050 

H^s 0.0278 ± 0.0041 0.0354 ± 0.0061 0.0341 ± 0.0053 

H^^ - 0.045 ±0.027 

H^^^ - 0.036 ±0.024 

iJf^5 - -0.009 ±0.023 

H^Q - -0.012 ±0.017 

Hf^.j - 0.020 ±0.013 

iJf'g - 0.004 ±0.019 

Hfj^ - 0.106 ±0.016 0.095 ±0.014 

i7£4 - 0.044 ±0.013 0.049 ±0.012 

if^s - 0.0343 ± 0.0089 0.031 ± 0.0087 

Hf^Q - 0.009 ±0.010 0.0068 ± 0.0094 

H^j - 0.0161 ± 0.0078 0.0195 ± 0.0069 

iffjg - 0.021 ±0.018 0.016 ±0.013 

Cp,49A 1.013 ±0.016 1.017 ±0.016 1.017 ±0.016 

Cd,49A 1.005 ±0.015 1.007 ±0.015 1.009 ± 0.015 

Cpa9B 1.023 ±0.015 1.039 ±0.016 1.031 ± 0.015 

Cd,49S 1.010 ±0.015 1.011 ±0.015 1.014 ±0.015 

Cp,6i 1.017 ±0.019 1.026 ±0.019 1.028 ± 0.019 

Cofii 1.004 ±0.018 1.018 ±0.018 1.018 ± 0.018 

Cp,87 1.019 ±0.015 1.029 ±0.016 1.025 ± 0.015 

Cd,87 1.008 ±0.015 1.007 ±0.015 1.011 ± 0.015 

Cp's9A 1.029 ±0.019 1.057 ±0.023 1.041 ± 0.020 

^D,89A 1.005 ±0.019 1.011 ±0.020 1.011 ± 0.020 

Cp,89B 1.016 ±0.015 1.024 ±0.016 1.021 ± 0.015 

^D,89B 1.003 ±0.015 1.004 ±0.015 1.007 ±0.015 

Cd'i39 1.010 ±0.015 1.012 ±0.015 1.014 ± 0.015 

X^/NDP 1274.3/1223 1248.0/1243 1255.6/1243 



Table 4: Correlation coefficients for the parameters of the ffi motivated by -R1990 parametriza- 
tion. 

61 62 h 

hi 1. -0.61 0.38 

62 -0.61 1. -0.87 

63 0.38 -0.87 1. 
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Figure 2: The values of H^ix) from the fits with model independent form of HT contribution 
to -Fi (empty circles - unconstrainted ffi, full circles - the fit with the constraint H^{x) = 
H[{x) For the better view the points are shifted to left/right along x-axis. The curves are 
predictions of IRR model for Q^ = 2 GeV^. 
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Figure 3: The same as Fig. 2 for Hl{x). 



We performed the fit using the QCD expression for F^ with the account of TMC 

2 T-,{P,D),LTf^ ^N 



■(6-2r) /" dz 



g2 ^2 V ^J^ ^2 

and additive form of HT contribution to F^, i.e. 

Fr^^^^ix, Q) = Ff '-)'™^(x, Q) + Hr\x)i^, 

where Hj^ ' (a;) is parametrized in the model independent form analogously to H2 ' (x) 
and /i2 ' (^)- The results of this fit are presented in column 2 of Table III and in Figs. 2, 3. 
One can see that with model independent parametrization of HT contribution to F^ the 
renormalization factors for old SLAG data noticeably increased. We remind in this connec- 



tion that in the source paper ||T2[ these data were renormalized using the data of dedicated 
SLAC-E-140 experiment. As far the latest did not reported data on proton target, renormal- 
ization of proton data was performed using bridging through SLAC-E-49B data. There is 
no possibility to conclusively choose between our renormalization scheme and the one used 
in [0. More proton data are necessary to adjust old SLAG results. As to deuteron data, 
their normalization factors practically do not deviate from 1. The errors of H^{x), due to 
SLAC-E-140 deuteron data, are significantly smaller, than for H[{x). In view of large errors 
of the last, HT contributions to Fl for proton and deuteron are compatible within the errors 
and we performed one more fit imposing the constraint H^{x) = H^{x). The results of this 
fit are presented in column 3 of Table IV and in Figs. 2,3. One can see that x^ obtained in 
this fit is practically equal to the value of x^ obtained in the fit without constraints minus 
the number of additional parameters. Comparison with Fig.l shows, that if Fl is fitted, the 
HT contribution to F2 slightly growth.. The value of strong coupling constant is insensitive 
to the constraint imposing. Their value 

as{Mz) = 0.1170 ± 0.0021(stat + syst), 

correspond to 

A^ = 337 ± 29{stat + syst) MeV, 

obtained at Q^ = 2 GeV^ with the help of relation 



as{Q) ^"^ 



2-K ln(21n(Q/A) 



/3oln(Q/A)L /5o/3 ln(Q/A) 
where 

and rif = 3. 

We checked how are the analysed data sensitive to the twist-6 contribution to Fl. For 
this purpose we fitted data with Fl{x, Q) expressed as 
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Figure 4: The dependence of twist-6 contribution to Fl on x. 



where functions Hj^ (x) and Hl{x) were the same for proton and deuteron and parametrized 
in model independent way. The resulting behaviour oiHf^ (x) is presented in Fig. 4. One can 
observe the trend to the negative values at highest x, but with poor statistical significance. 
The x^ decrease in this fit is about 25 (remind that standard deviation of x^ is V2 ■ NDF, 
i.e. is about 40 in our analysis). These results are in correspondence with the estimates 
of twist-6 contribution to R{x, Q) presented above - the fitted twist-6 contribution to Fl 
is at the level of one standard deviation. In our studies we did not completely accounted 
TMC correction of the order of 0{M^/Q^) and then, in a rigorous treatment, Hj^ does not 
correspond to a pure dynamical twist-6 effects. Meanwhile the contribution of the omitted 



r(4) 



terms to if| estimated using formula from |]T5[, do not exceed 0.001 in our region of x and 



hence can be neglected. 

3. For the comparison with other parametrizations we calculated deuteron R{x, Q) using 
the relation 



R 



D,HT, 



X,Q) 



Fr''{x,Q) 



FF'^'^ix^Q) 



:^D,HT, 



X,Q) 



1 + 



4M^x 



2^2 



Q' 



and the parameters values from column 3 of Table III. The obtained values of R are presented 
in Fig. 



5 together with -R1990 [H 



and -Ri998 113 parametrizations. In average all three 
parametrizations coincide within errors, although our curves lay higher at the edges of x- 
region and lower in the middle. 

The same tendency is valid for HT+TMC contribution to R, evaluated as the difference 
between R^'^'^{x,Q) and R^'^^{x,Q), where 



:^D,LT, 



R 



D,LT, 



X,Q) 



x,Q) 



:^D,LT, 



X,Q) 



Ff'^^(x,Q)' 



This contribution is presented in Fig. 6 together with the power parts of the variant b) of -R1990 
and -Ri998 parametrizations. Since error bands are not available for latest ones, only central 
values are pictured. Relative difference between log parts of -R1990 and -R1998 parametrizations 
and R^'^^{x, Q) is very large at highest x, although absolute difference is not significant due 
to smallness of these terms. In Fig. 8 we present the TMC contribution to R, calculated as 
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Figure 5: One standard deviation bands for our value of R^'^^{x,Q) (full lines), i^iggo 
(dashed lines) and -Riggs (dotted lines) at Q^ = 2 GeV^ Our value of R is obtained as spline 
interpolation between points x = 0.3, 0.4, 0.5, 0.6, 0.7, 0.8. 
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Figure 6: One standard deviation bands for HT+TMC contribution to R (full lines). Dashed 
and dotted lines correspond to power contributions to i?iggo and i^iggs- Dashed-dotted line 
represents the power term of BRY parametrization [14]. 
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Figure 7: One standard deviation bands for LT contribution to R (full lines). Dashed and 
dotted lines correspond to log contributions to i^iggo and -Riggs- 
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Figure 8: One standard deviation bands for HT contribution (full lines) and TMC contribu- 
tion (dashed lines) to R. Lower and upper bands for TMC are practically indistinguishable. 
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Figure 9: The E140X data on nucleoli R{x,Q) (empty circles). Inner bars correspond to 
statistic errors, total bars - statistic and systematic combined in quadratures. Also are 
presented calculations of deuteron R{x, Q) performed on our fit (full circles). For the better 
view the points are shifted to left/right along x-axis. 



Table 5: Coefficients of correlations between the HT parameters and as{Mz)- Figures in 
parenthesis are global correlation coefficients for the HT parameters. 



X 

0.3 

0.4 
0.5 
0.6 
0.7 
0.8 



H^[ 



X] 



H. 



D, 



X] 



HAx) 



-0.507(0.914) -0.707(0.964) -0.029(0.868) 

-0.847(0.955) -0.910(0.976) 0.122(0.852) 

-0.909(0.968) -0.944(0.987) 0.244(0.870) 

-0.905(0.971) -0.917(0.977) 0.222(0.866) 

-0.871(0.972) -0.901(0.980) 0.072(0.894) 

-0.460(0.868) -0.429(0.875) 0.156(0.870) 



^D,TMC _ j^D,LT ^j^g^g 



R 



D,TMC, 



^D,TMC, 



X,Q) 



X,Q) 



;,D,TMC. 



x,Q)-Fl 



D,TMC, 



X,Q) 



1 + 



AM^x'^ 



Q' 



1, 



and the HT contribution, calculated as R^'^^ — R^''^'^^'-^ . It is seen that the TMC contribu- 
tion also is significantly smaller at large x than the HT one, so the latest is dominating in 
this region. Meanwhile we should note in this connection that as it was observed in p, p!7| , 
account of NNLO QCD can diminish the value of HT contribution. This effect can be seen 



in Fig. 6, where we present power part of R parametrization ||T^, obtained with the partial 
account of NNLO. 

In Fig. 9 data of SLAC-E140X experiment on nucleon R{x, Q) [|18|, which were not included 
in the fit, are compared with our R^{x, Q), calculated at the parameters values from column 
3 of Table HI. One can observe a good agreement of the data and our parametrization. 
We also compared the results of model independent analysis with a predictions of infrared 
renormalon model (IRR) [^ . In this model the HT x-dependence is connected with the LT 
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x-dependence. In particular for nonsinglet case twist-4 contribution is expressed as 

Jx 

C2iz) = - ^ +2{2 + z + Qz^) - 95(1 -z)- 5'{l - z), 
Cl{z) = 8^2- 45(1 -z), 

Po 
where Cp = 4/3, C = —5/3. The normahzation factor A/j can be considered as a fitted 
parameter, or, in other approach, set equal to the value of Aqcd- In Figs. 2,3 we present the 
IRR model predictions for Q^ = 2 GeV^. At this value of Q the number of active fermions 
Uf = 3 and we set Ar = A-^ = 337 MeV, as obtained in our analysis. Parameters of 
LT structure functions were taken as in column 2 of Table III. One can see that the model 
describes the data on Hl{x) fairly well. This in agreement with the curves, presented in 
Fig. 6 of Ref. [|y] (remind that our value of A is about 1.3 times larger, than used in |2^ 



to calculate these curves). At the same time there is evident discrepancy between the model 
and the data on H2{x). In this connection we should note that H2{x) is strongly correlated 
with other parameters, and in particular with as{Mz) ■ In the Table V we present a global 
correlation coefficients pm, calculated as 



Pm y -L -L / (^mm / ^n 



where Cmn is parameters error matrix, e^n - its inverse and indices m, n runs through fitted 
parameters. Global coefficients characterize the extent of correlation of a parameter with 
all others parameters. (In the case of the fit with two parameters pi = p2 and equal to the 
correlation coefficient between these parameters). From the Table V one can see that Hi 
is almost uncorrelated with as{Mz) and is correlated with other parameters less than H2. 
This fact can be readily understood qualitatively. The total value of R^^ is defined from 
the y-dependence of cross sections and is weakly correlated with as- Then, the correlation 
of HT contribution to Fl can arise only due to dependence of LT contribution on as- As 
far, the LT contribution is small comparing with the HT one (see Figs. 7,8) this dependence 
does not cause significant correlations of the HT contribution and as- Due to lower values 
of p, Hl is more stable in respect with the change of the ansatz and input of the fit. The 
change of as towards higher values would not affect Hl, but can decrease H2- 

Conclusion 

In conclusion, we performed NLO QCD fit of the combined SLAC/BCDMS/NMC DIS 
data at high x- Model independent x-shape of high twist contribution to structure function 
Fl is extracted. Twist-4 contribution is found to be in fair agreement with the predictions of 
infrared renormalon model. Twist-6 contribution exhibits weak tendency to negative values, 
although in the whole is compatible with zero within errors. 
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